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ABSTRACT 


This  is  the  final  report  on  investigations  of 
theoretical  problems  in  cosmical  gas  dynamics.  An 
attempt  has  been  made  to  derive  a  kinematical  model 
of  the  Orion  Nebula  from  observational  models  of  the 
density  distribution  and  data  on  the  internal  kine¬ 
matics.  Numerical  methods  have  been  developed  for 
integrating  the  equations  governing  the  dynamics  of 
H  II  regions  by  the  method  of  characteristics. 
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THE  EXPANSION  OF  THE  ORION  NEBULAE 


This  investigation  was  an  attempt  to  construct  a 
kinematical  model  of  the  Orion  Nebula  on  the  basis  of 
1)  Spherically  symmetric  models  of  the  distribution  of 
the  electron  density  derived  from  optical  and  radio 
observations  (Osterbrock  and  Flather  1959;  Menon  1961) 
and  2)  observational  data  on  the  internal  kinematics  of 
the  nebula  (Wilson,  I-’unch,  Flather,  and  Coffeen  1959). 
The  nebula  was  treated  as  an  expanding,  spherically- 
symmetric  system  in  which  random  motions  of  the  gas  are 
superimposed  on  the  velocity  of  expansion.  Under  such 
conditions,  the  average  radical  velocity  <V(s)> 
of  the  gas  observed  in  a  given  spectral  line  along  a 
line  .of  sight  at  a  projected  distance  S  from  the  center 
^  of  the  nebulae  is  related  to  the  velocity  of  expansion 
UU|  ,  where  A  is  the  true  distance  from  the  center, 
in  the  manner 


JsCh)  a  cn.)  esi  9  e  <Ll! 


<V(s)>  = 


/6(A)  6~r  <U 


(1) 


Here  the  integrations  are  performed  along  the  line  of 
sight  through  the  nebulae , is  the  optical  depth  into  the 
nebula,  6(a)  is  the  volume  emissivity  in  the  spectral 
line,  and  the  factor  cos <9  projects  the  velocity  of  expansion 
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onto  the  line  of  sight.  Values  of 

were  derived  from  the  catalogue  of  Wilson,  et  al. ,  while 
values  of  6(A)  were  derived  from  the  papers  of  Osterbrock 
anc*.  Flalher  and  of  Kenon.  Optical  depths  into  the  nebula 
were  derived  from  results  of  an  investigation  by  O'Dell 
and  Hubbard  (19! 5,  private  communication).  The  objective 
was  then  to  find  plausible  functions  U  (A)  such  that  equa¬ 
tion  (1)  would  be  satisfied  over  the  range  of  $  for  which 
<V(5)>  could  be  determined. 

For  the  model  envisaged,  the  expectation  was  that  t-tlA) 
should  vanish  at  A  =  0  ,  increase  to  a  maximum  value  IKo 
(say)  at  some  distance  -  A,  (say) ,  and  fall  to  zero  as 
00.  Dynamical  considerations  imply  that  llo  .’>ould  be 
of  the  order  of  the  isothermal  speed  of  sound  (about  10 
km-sec"1)at  most.  It  was  found  that  equation  (1)  could  be 
x  satisfied  only  if  Ot(A)  attains  a  maximum  value  of  the 

order  of  several  hundred  kilometers  per  second.  This  result 
is  unacceptable  because  it  violates  the  aforementioned  dyna¬ 
mical  considerations  and  because  it  implies  that  the  dispersion 
of  the  radial  velocities  in  the  nebula  should  be  much  larger 
than  is  observed. 

When  this  result  was  obtained,  it  was  recognized  that  one 
possible  source  of  error  in  the  analysis  lay  in  the  assumption 
that  the  Orion  Nebula  cluster  is  at  rest  relative  to  the  center 
of  the  Nebula.  On  this  basis,  ^  V  (S)^  had  been  determined 
relative  to  the  cluster.  This  procedure  would  give  erroneous 


results  for  U(A)  if  the  nebula  were  drifting  relative  tc 
the  cluster.  Subsequently  studies  of  the  radio  recombina¬ 
tion  lines  of  hydrogen  revealed  such  a  drift  (Mezger  and 
Hoglund  1967) .  Departures  from  spherical  symmetry  must  be 
additional  sources  of  error.  Indeed,  more  recent  investi¬ 
gations  (e.g.  Webster  and  Altenhoff  1970)  nov?  indicate  that 
the  Orion  Nebula  is  much  less  regular  and  symmetric  than  was 
believed  when  this  work  was  undertaken.  For  these  reasons, 
it  now  appears  that  a  reliable  kinematical  model  cannot  be 
derived  along  these  lines. 


NVr-.IRICAL  METHODS  FOR  THE  INVESTIGATION 


OF  THE  EVOLUTION  OF  II  II  REGIONS 

It  had  been  proposed  under  this  Grant  to  investigate 
aspects  of  the  early  evolution  of  an  H  II  region  which  had 
become  of  interest  in  earlier  work  (Vandervoort  1966) .  In 
pa  "ticel-ir ,  it  was  intended  to  consider  further  the  effects 
of  f  n  .  e  time  scale  for  the  formation  of  the  exciting 
star  to  include  a  more  careful  treatment  of  the  ioniza¬ 
tion  and  heating  of  the  g  is  in  the  formulation  of  the 
dynamics  of  the  H  II  region.  A  central  question  was  whether 
the  early  evolution  of  an  H  II  region  would  be  characterized 
by  an  ionization  front  of  the  R-type  propagating  directly 
into  the  undisturbed  neutral  gas  or  a  front  of  the  D-type 
following  a  shock  into  the  undisturbed  gas.  And  if  the  front 
were  of  the  D-type,  the  question  arose  whether  or  not  a  rapid 
increase  of  the  luminosity  of  the  exciting  star  would  induce 
a  transition  to  a  front  of  the  R-type.  This  question  had  been 
raised  earlier  in  an  investigation  of  non-stationary  ioniza¬ 
tion  fronts  (Vandervoort  1965  a,b) . 

It  was  clear  at  an  early  stage  of  the  work  that  ';he 
development  and  use  of  numerical  methods  would  be  required. 
Accordingly,  the  thrust  of  the  work  shifted  to  the  development 
of  suitable  numerical  schemes  for  integrating  the  basic  equa¬ 
tions.  At  the  time  this  work  was  undertaken,  other  investigators, 


v  •>o^tI^«S*iS%5Ss5^!!*e=J',:B: 


particularly  Mathews  (1965)  and  Lasker  (1966)  ,  had  investi¬ 
gated  the  evolution  of  11  II  regions  by  numerical  integration 
of  the  hydrodynamic  equations  in  a  Lagrangian  representation. 

It  was  thought  that  the  present  work  might  compliment  those 
investigations  if  it  were  based  on  other  numerical  schemes. 

For  this  reason  and  because  the  Principal  Investigator  had 

had  previous  experience  in  the  use  of  the  method  of  charac- 

\ 

teristics  (Courant  and  Friedrichs  1949) ,  it  was  decided  to 
develop  numerical  methods  for  integrating  the  basic  equations 
by  the  method  of  characteristics. 

In  particular,  such  a  scheme  has  been  developed  for  the 
numerical  integration  of  the  equations  governing  the  structure 
of  a  non-stationary  ionization  front  (Vandervoort  1965a;  all 
equation  numbers  in  this  section  of  this  report  refer  to 
equations  in  that  paper) .  Characteristic  equations  corres¬ 
ponding  to  equations  (1),  (2)  ,  (4),  and  (6)  we re  derived  in 
a  conventional  manner  (Courant  and  Friedrichs  1949) .  These 
were  replaced  with  a  system  of  finite  difference  equations 
appropriate  for  a  numerical  integration  by  a  "method  of 
specified  time  intervals"  (Lister  1960) .  At  each  time  step 
in  the  integration,  the  flux  must  be  determined  by  integrating 
equation  (3) ;  this  was  accomplished  with  the  aid  of  a  conven¬ 
tional  scheme  for  numerical  quadrature.  The  boundary  conditions 
expressed  by  equations  (8)  -  (10)  were  replaced  with  boundary 
conditions  at  the  first  and  last  points  on  the  spatial  grid. 
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A  program  for  the  execution  of  this  method  has  been 
written,  and  it  has  been  tested  on  the  problem  of  a 
stationary  ionization  front  for  which  an  exact  solution 
is  known  (see  §IV  in  Vandervoort  1965a) .  Considerable 
experimentation  was  required  to  find  satisfactory  combina¬ 
tions  of  the  number  of  points  in  the  spatial  grid  and  the 
size  of  the  time  step.  The  status  of  the  work  at  the  end 
of  the  grant  period  is  that  it  has  been  possible  to  pro¬ 
duce  accurate  solutions  for  the  structures  of  stationary 
R-type  fronts  for  intervals  of  time  greater  than  the 
character' stic  time  scales  during  which  elements  of  the 
neutral  gas  are  over-run  and  completely  ionized.  These 
results  encourage  the  belief  that  with  further  development 
these  methods  could  be  applied  to  non-stationary  fronts  and 
generalized  for  application  tc  other  aspects  of  the  dyna¬ 
mics  of  H  II  regions.  In  view  of  the  circumstance  that 
such  work  would  involve  a  considerable  investment  of 
computer  time,  for  which  funds  are  currently  unavailable, 
work  on  this  project  was  suspended  at  the  end  of  the  grant 
period. 
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